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ABSTRACT 


Decreasing price-performance ratios have resulted 
in a large number of minicomputer installations in the process 
industries. This thesis describes multiloop control software 
for a Data General Corporation Nova Computer. 

The software is developed for digital (D(z)) and 
proportional-integral-derivative (PID) controllers intended 
for use in undergraduate education. The machine language 
encoding of the program allows the user to study sampled-data 
feedback control systems in a convenient manner. 

Hybrid computer simulation of several digital control 


systems is presented to illustrate the use of the system. 


Pa 


aeacoxrg Bie ak coat sAtal Bs 
arawli0e Loxsn02 qooLtstum’ sod 
bas ((3)@) testa =e eqoL: 


babassat nselonaigs (a) 2 


Sy 


ACKNOWLEDGMENTS 


The author wishes to thank Professor Y.J. Kingma, 
thesis supervisor, for his encouragement and guidance in the 
preparation of this thesis. The author also wishes to thank 
other staff members and graduate students in the Department 


of Electrical Engineering for their assistance. 


The financial assistance provided by the National 
Research Council and Department of Electrical Engineering is 


gratefully acknowledged. 


oy uh a lt a i vi 
- asi3 ak sonsbivg bow 40 
th taal 


Jnedt 03 asdetw oals 


| angmtzeqad sf -ak ean 


CHAPTER 1 


CHAPTER 2 


CHAPTER 3 


CHAPTER 4 


REFERENCES 
APPENDIX A 


APPENDIX B 


TABLE OF CONTENTS 


INTRODUCTION 
CONTROL SYSTEM IMPLEMENTATION 


2.1 Organization 

2.2 Interrupt Structure 

2.3 General Program Routines 
2.4 D(z) Implementation | 
2.5 PID Implementation 


SYSTEM TEST RESULTS 


3.1 D(z) Controllers with Common Sampling Rates 


3.2 D(z) Controllers with Individual Sampling 
Rates 


3.3 PID Controller 
3.4 PID Controllers with Common Sampling Rates 


3.5 PID Controllers with Individual Sampling 
Rates 


3.6 Calculator Mode 
USER EXPERIENCE AND CONCLUSIONS 


4.1 User Experience 


4.2 Conclusions 


20 
20 
28 


28 
34 


34 
37 
40 


40 
42 


43 
44 


61 


/ 


ested golignst nome? dyiw ezelle 
gniiqus® leubtvibal dite nn were ml = 
reflow) Gs Ge 
agiad goblqase gomeod tbe enOrLgAAED, GS Ae. 
galiquse fsubtytbal Adiw asefforsi0') ara Z.€ 
soten 


2% £88 2 a 


- sbol rossSieotad 8.€ 
eMOTeUZIMOD CMA tase aueu 


Os 
OA 
Ss 
eA 
ba 
fd 


LIST OF TABLES 


TABLES Page 


321 Tabulation of D(z) Controllers Used in 
Figures 3.5 and 3.6 29 


4.1 Loop Execution Times 41 


LIST OF FIGURES 


FIGURES Page 
bel Computer Control System 3 
ae System Configuration 4 
Zak System Organization 7 
Zee Program Interrupt 10 
2e3 Timing Sequence for Control Loops Lt 
one Digital Controller 16 
225 Realization of D(z) Transfer Function 16 
2-6 PID Control System 18 
vb Digital PID Controller 19 
sega Response of Systems with Second Order D(z) 

Controllers 21 
3.2 Controller Response for Systems in Figure 3.1 M9) 
SE Response of Four Different Systems 23 
3.4 Controller Response for Systems in Figure 3.2 24 
335 Response of Systems with 2,3,4 and 5th Order D(z) 

Controllers 26 
3.6 Controller Response for Systems in Figure 3.5 27 
Sar Response of Systems with Sampling Rates of 0.1, 

0.4, 0.7, and 2 Seconds 30 
Shs Controller Response for Systems in Figure 3.7 31 
339 Response of PID Control Elements to a Step Input 32 
3.10 Response of PID Control Elements for the System 

1/(2s+1) (10s+1) 33 
3eh1 Output Response of 1/(2s+1)(10s+1) for Various 

Control Settings 35 
Sal2 Output Response of 1/(2s+1)(10s+1) for Various 


Sampling Rates with Control Parameters Fixed 36 


if 


of 


meteye lorinod xSauqitoD — 


nobismsiprod RGIS 


Bgooi fox41n09 10% Somgupse or: 
xefforine® lettgid 
aotjaxtiaall 
fousacd: arg 
ot keraasD aQqq Inaigtd 


aneteve’ to.eanogean 


(s}d x9bx0 baosee dtiw 
erelloyinod 


I.€ s1gtl mt emsjaye x02 senoqasd +alforjae) — 
amas2ye Jnstetiid r0F to senogeal \ 
S.€ stugtt ak one ER2 x02, senegaal weligxind 


Cs ie tob10 d3@ bas 4,€,8 daitw ains3eye ‘to s senoqaga 
 atsLiortnod 


Gat suugit al amjeye 192 seaoga’i soltorson 


,1,.0 to Loans Bake qm & Ast emoaaye'to' 6 
' ebnoasé s bas ee nd 


T.€ stugh® nt onaraa ret aenoqaoh eae F a 


tuqnel asi 2 03 enema la tox3¢00) asa bi hint 


S+a0L) (i-+e8) opr 
ba wscer ditwie 


Response of Four Different Systems and Sampling 
Rates of 0.5, 0.1, 1.0, and 2.0 Seconds 


Calculations in Calculator Mode 

Generate Operating System 

Restart Program 

Interrupt Handler 

D(z) Control Loops with Common Sampling Rates 
D(z) Control Loops with Individual Sampling Rates 


Calculates D(z) Values to be used for the Next 
Sampling Period 


Common Sampling Rate PID Loop 
Individual Sampling Rate PID Controller 
Multiply Routine 

Overflow Check Routine 

Message Type Routine 

Decimal to Octal Conversion 

Floating Point Input 

Calculator Mode 

Parameter Storage Tables 

Start-up Operating Message 


Hybrid Computer System 


38 
39 
46 
47 
48 
49 


50 


a1 
By 
53 
54 
oy) 
56 
57 
58 
59 
60 
87 


88 


vB 


88 


-_ 
y : uv 5 
=n nm! j 7 al 
’ is * ae! , 7 
4 4 7 : au \ r ai ies be 
F ‘ , ¢ ‘ t a “s Ni i 
j n ‘ : 4 say iT re 
is , ; . = It's 
i [eo A oy 1 7 7 


aattamed bas aingd2¥2,ans78td Li awoil Jo sehoqeeh a a 
Bhnonee wy qt na 2D. t am oe 0 +e ) te eenaf Ui, 
pi > a 


aban toda Led nt énbiitetbogan 
5 ee 


nstaye: SnitaerTsqo. sisrenad 


“I 5 Ol (“7 a 
fiasgett jasgesh 8 | Rell 


toibash iquirwsant oa 
t 


ests gnifquse romnod: dtiw Boos Latzae) (s)d 


asse1 antiamse Laub ibal datw aqoad igi (s)a = Cel 


7x9 si3 yo? bsav sd od eentey ts)a gotiebuo tat 4 
hobte? galfqmie ; 


“a! 


gous dit sia# BR) tL qaas gomusoD 
is{foraa09 GId saan sek toaie LaubivPbal 
spiguod yiqkiivM a A 


gmtivoy slostid woltsevd . 
antiuad 'sqyt sg siete 


noletavaod feand of thetasa al 


Sarqerd takes anne 


abo sorssied od 
r ” i it ie We 
eatdnT susie ail id ae 


ae 
egsedelt sph eae quae ed 
An 


| tan 


ie 
ngjaye r93uqmod bird 


q 
- ni 
Re. 3 a Ww 
’ ; | : 


CHAPTER 1 


INTRODUCTION 


During the past several years digital computer technology 
has experienced rapid growth. The small computer is now becoming a 
useful and economical device in such areas as system simulation, on- 
line data processing, and on-line control. The undergraduate engineer 
working in the systems area must be capable of working with digital 
system concepts. 

The introduction of general purpose computers has made 
possible real-time digital control systems for industrial applications. 
There are two approaches available for introducing digital computers 
into undergraduate control education. They are: (1) require the 
student to write his own programs for solving problems or (2) allow 
the student to use existing programs to solve specific problems. The 
first approach requires large amounts of time and programming effort 
and the second method, while not allowing him to follow each step, 
permits the student to study the effectiveness of a digital computer 
solution to the problem being studied. 
fibevenw! Mfg shown that a digital control laboratory for 
undergraduate instruction is a desirable feature of the undergraduate 
ciriculum. The student becomes familiar with discrete time events, 


and develops an appreciation of the interaction between the digital 


computer and external devices. The student is familiar only with the 


* Numbers in [ ] refer to references listed at the end of the thesis. 
References are listed in the order in which they appear in the text. 
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digital computer as a general-purpose calculator from his preceding 
undergraduate courses. However, using the digital computer in a 

real time control application, will acquaint him with it as an elect- 
rical device, using signals on a binary level, analog level or in 
terms of a pulses. Martens has found that the undergraduate is 
capable of understanding digital concepts and that the graduating 
engineer will be familiar with the tasks involved in development of 
computer control systems. 

The purpose of this thesis is to develop an operating 
system for a Datagen Nova computer to be used in an undergraduate 
control laboratory. One digital computer is used in a sequential 
time-sharing mode for all control loops. The advantages are that this 
permits a central supervisory control scheme, multivariable control, 
and simultaneous operation of several interdependent loops. Since 
the control laws are implemented through software, unlimited flex- 
ibility exists in the implementation of the control functions. 

The Wova used for this system has 8K of memory, a real- 
time clock, hardware multiply/divide and multichannel Analog to 
Digital (A/D) and Digital to Analog (D/A) converters. Although the 
Nova software contains a BASIC compiler and a Floating Point instruc- 
tion set, they are of limited usefulness for laboratory studies. It 
is therefore necessary to restrict all programming to the computer's 
symbolic assembler language. This language provides all the cap- 
abilities needed in the laboratory. 

The operating system consists of sampled-data control func- 


tions of two types. These are the z-transform (D(z)) controller and 
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and a digital Proportional-Integral-Derivative (PID) controller. The 


control system is shown in Figure 1.1. 


e, (t) : controlled 
Hold output 


reference 


Digital 


r(t) Computer CircuLt 


Figure 1.1 Computer Control System 


The operating system is generated via the teletypewriter and can be 
tailored to the number of analog computers on-line and to the experi- 
ment being studied. The program is self-starting, thus allowing 
operation by those who are not necessarily familiar with machine 
language programming. The functioning system allows the perdent to 
operate individual analog computers simultaneously without regard to 
other operations being performed at that time. The teletype handles 
their requests for changes in control constants. The input of these 
constants (the foreground program) does not disturb the operating 
control program (the background program). The system configuration 
is shown in Figure 1.2. 

Some of the specific design requirements for the program 


were that it should: 
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Figure 1.2 System Configuration 


a. Be simple and convenient to use but also be capable of 
handling realistic problems. 


b. Handle standard control functions as selected by the user. 


c. Allow easy expansion of the program to service additional 
on-line devices. 


d. Be capable of sampling rates of .1 seconds for ten processes. 


e. Generate an output in minimum time after an input sample is 
applied to the control loop. 


f. Run in a time-sharing mode on the Nova computer. 


g. Occupy a minimum of core memory. 


The final program is described in Chapter 2. A general 
description of the program is given by Figure 2.1. Flow charts for 
the main subroutines are given in Appendix A and the program listing 


with comments is given in Appendix B. The methods used in formulating 
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a fast-multiply routine, modification of floating point numbers, and 
calculation routines are also discussed in Chapter 2. The methods 
of realizing D(z) controllers are given in section 2.4. The PID 
implementation and methods of integral and derivative control are 
given in section 2.6. 

Chapter 3 is concerned with the operation of the system. 
Test results are given for the simultaneous operation of four analog 
computers. Also shown is the response of the PID elements in execu- 
tion of a control function. 

Chapter 4 includes conclusions and a discussion of user 


experience. 
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CHAPTER 2 


Lea Organization 


The basic system organization is shown in Figure 2.1 Upon 
entry of the program, storage areas for controller parameters, delay 
elements, timing tables, mode tables, and calculation areas are 
initialized to avoid problems arising from random memory words in those 
locations from a previous user's program. After completion of this 
function, the program requests the number of analog computers or 
processes on-line and then requests a number representing a desired 
system operation program. Upon specification of the type of program 
desired, the program generates the necessary links to the real-time 
clock and teletype for that operating system. After completing this 
task, the interrupt is enabled and the system is operating in the 
specified mode. At this time the user can define control parameters, 
the sampling rate, or enter the calculator routine. 

The system can only be changed by manually stopping the 
computer and restarting the program in memory location 000002. This 
allows "hands off" operation of the computer in that all requests are 
handled via the teletype and the computer will not allow invalid 
requests which could disrupt the running system. 

Control parameters are entered in floating point format and 
are converted to a modified flcatins point representation thus allowing 
faster multiplication by dealing with numbers in a smaller magnitude 


range. The sampling rates specified by the user must be multiples of 


the clock frequency. 
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A calculator option also exists which permits the basic 
mathematical operations of addition, subtraction, multiplication and 
division. This routine can be entered at any time since it is run 


as a foreground program. 


2.2 Interrupt Structure 

The running program is operated entirely by the interrupt 
facility on the Wova computer. The input-output devices must be 
serviced infrequently relative to the processor speed and only a small 
amount of processor time is required to service them. These devices 
must be serviced within a short time after they request it. The 
program interrupt allows interruptions in the current program sequence 
and allows concurrent operations of the main program and a number of 
peripheral devices. The hardware allows conditions internal to the 
processor to signal the program by requesting an interrupt. The 
interrupt requests by a device are controlled by its Done and Interrupt 
Disable flags. A device completing an operation sets Done and this 
action initiates a program interrupt if Interrupt Disable is clear. If 
Interrupt Disable has been set by the program, the device cannot request 
an interrupt. 

The processor starts an interrupt after completion of the 
current instruction, with interrupts enabled and a device waiting for 
an interrupt to start. The processor clears Interrupt On, saves the 
Program Counter in location 0 and simulates a JMP @ 1 to enter the 
interrupt service routine. 

The interrupt service routine determines which device requires 


service by testing with I0 skips. The program servicing the device can 
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enable interrupts and establish a priority structure that allows higher 
priority devices to interrupt the current device service routine. In 
this case, location 0 must be saved along with the carry and those 
accumulators which are used in the device service routine. Priority 

is determined by a mask that controls the states of the Interrupt Dis- 
able flag in the various devices. Any device whose Disable flags is 
set cannot cause an interrupt by setting of Done. Those devices in 
which the Interrupt Disable is clear can interrupt the current routine 
and are regarded by the processor as being of higher priority. 

After completion of a device service routine, the routine 
restores the pre-interrupt states of the accumulators and carry, turns 
on the interrupt and returns to the interrupted program. If the ser- 
vice routine allows interrupts by higher priority devices, then before 
dismissing the current routine, the routine will turn the interrupt 
off to prevent further interrupts during dismissal. In dismissing, the 
routine will re-enable lower priority devices that were not allowed to 
interrupt the current routine but will be allowed to interrupt the 
program to which the processor is returning. The processor always 
executes one more instruction (return to the interrupted program) before 
another interrupt can start. 

The program listed in Appendix B uses the interrupt facility. 
In the case of deadbeat sampled-data controllers, it is necessary to 
perform operations on time. The real-time clock is made the highest 
priority device and it controls the timing of all events occurring in 
the background program. The teletype has a priority level of one below 


that of the clock, thus enabling it to handle all input/output requests 
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at the fastest possible rate. Figure 2.2 indicates a typical sequence 


of events during program operation. 


2.3 General Program Routines 
2.3.1 Program Sequence 
The real-time clock controls the sequence of events result- 
ing in control outputs. The case where all control loops are operat- 


ing with a common sampling rate is shown in Figure 2.3. 


t : t 
Service Calculations Service 
Required Before Next Required 
Sample 
Control Control Control Control 
0 n 0 n 


T = period between samples 


Figure 2.3 Timing Sequence for Control Loops 


Each control loop will have a fixed sampling rate and a 
constant period T between samples. When the system is operating with 
each loop having its own sampling rate, software clocks are provided 
which time-out at the required intervals as specified by the user. 


In this case, idle loops are entered when no output is required and 


a : aan ts 


somsupse Isoltgy) s setaokbut S.S syygtl ater efdvesoq. saaxeat of 
: ; } » 2 a 
. «tO fisitsgo metgotq gainub = TSN 


aap tuo | maz 272 cera. es 
| | ce 


Sareyppes mstgord ha c.8 \ 


see 
; if 
-glueex etasve to soneupee ada efowsmoo doef smit~ises sat - me 


“Pa 


‘—tayeqo ots aqool foxnoo Iie stsdw sens ent , .2eduqdun Lowines | 


-£.8 aaugt T ot mwosde el stv gaklqmss commoo # dt bw 
4 
} v 
$+ ey a mene 
i ree } ee | 
Ran nn ee te me = 
t eho) : i BLA 4 4\ i] 
-~ a — Te wy — emer ry -———— a. oe v 
| | | | 7 
i j } ' ' 
#; ' 
: ‘ , 
j ; ) ‘ 4 
: ; i ' 
‘ a ; 
al ate aneueenteed Hf 
ssivise enolisiuoits® | | ' 
bs1tupext tus Vi stole . 7 
glgqmsey ) oat 
Wn . 
fossq09 Lortao0D \ ). forded — iy 
fr 0 e ; | aT 4: wT 
¥ 7 il Fe 7 rr ie 
ae fq jms noo udurt bo.tas fa? 
7) a ¢ 
a ~, ; a A i “ols 
| wil ee, eee 
bs 7 hare * Lae a) : ¢ 
| aqoo! Loxsnep ‘toi fore ania £8 paged a : 
7 ao " ny wee | 7. : 8 7 ; 
an iY : ; ae | : ‘ i 4 
pe ee, & 
7 8 brs suey gat Lqmee be axtt 5 veil thw spot trom ta ty cae 
a ; | i 


ah _ 
Halts gatts st9q0 et bathe att naib spose 18% soso. poke aq 2 ae 309 


' (nat sire | * 


i | 

ye a tee 9x3 ailoots stew ioe oes gag 

ie stem, lt "ed" baiaioana eB 8 ‘slowsegnl he ht A 
"7 ao f * - ‘ i 


—- in - : i fine 7: . , og oy vies | 
| . _ betiupss auqibe ai new. birahad : ou 
4 . ; > ae 
ie F 
. ff 


12 


thus the period T for any loop remains constant regardless of the 
number of operations required after a clock interrupt. 

In the case where each loop has a different sampling rate, 
all loops must be complete before another clock interrupt can occur. 
This is not the case in the loops with a common sampling rate, since 
the clock can interrupt the control loop service program to increment 
a counter for lower sampling rates. This allows more devices to be 
placed on line as well as providing more time for other desired features. 
Appendix A illustrates the flow of the above two sequences. 

2.3.2 Modification of Floating Point Input 

A routine to modify the floating point representation of 
parameters is described next. Using the Wova floating point roariner 
the user can enter a value in engineering notation in the following 
form: 


Cem fs. nl* in ws ne tem im) pebreak: 


where characters in [ ] are optional entries. This number is stored 


in internal floating point as a 32-bit binary string of the form: 


where "S" is the sign of the mantissa "M''. The mantissa is considered 
to be a normalized six digit hexidecimal fraction with a range of 


magnitude of 


1 6 


16<M< 1-16 
The characteristic "C", is the integer exponent of 16 in excess 64,0 
code. This gives a range of magnitudes for floating point representa- 


tion of a decimal number as 
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8.7x 10/8 < ¥F ss ae Es 


Multiplication of two floating point numbers requires eight 
milliseconds which precludes using floating point instructions for the 
calculations needed in the control loops. Therefore a new representa- 
tion was formulated to handle numbers in a limited range. The mantissa 
was contained in a 16-bit word in the form of a four digit hexidecimal 
fraction. The characteristic and sign are contained in another 16-bit 


word. The useful range of this number system is 


which becomes 1.5 x liges SOPs Aa x ae This was found to meet all 


the requirements for the accuracy needed in a digital process control 
loop. Entry of a number less than 1.5 x 107 will result in a zero 


3 will result in it being 


answer and using a number larger than4.1x10 
represented by Hote 10 The routine which performs the above is 
called FLPT. 
2.3.3 Multiplication Routine 

In the calculation of output values, we deal with analog to 
digital converted values which are 12-bit binary numbers. Multiplying 
these values by floating point numbers requires a routine that will 
execute this function, check for overflows and, most important, complete 
this task in the minimum possible time. The analog-to-digital converter 


scales +10 Volts to 003777 and -10 Volts to 174000 with a resolution 


of + the Least Significant Bit (L.S.B. = 4.8 mV). 
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The multiply routine uses the hardware multiply device to 
multiply the signal value by the mantissa of the decimal number. 
Then the required number of shifts are performed as specified by the 
characteristic of the decimal number. The last operation performed 
is to affix the correct sign to the 16-bit result. This routine 
requires a maximum of 170 usec. The shortest time through this routine 
is 65 usec. 

2.3.4 Calculator Routine 

A calculator routine is provided, and it can operate during 
the idle time of the program. The routine provides for addition, sub- 
traction, multiplication and division (+, -, *, /). It is entered by 
typing CTRL/C and exited by typing the same command. An up arrow (‘t) 
indicates that the user is in this mode. 

Numbers are entered in engineering notation in the following 
sequence. 

x[breakl]y[break2] where breakl is the desired operation 
and break2 is a carriage return. The result is returned in engineer- 
ing notation. The range of numbers that this routine will handle is 


<7 
approximately 10 : 


2.4 D(z) Implementation 

There are four methods of realizing digital controllers: 
direct programming, cascade programming, parallel programming, and 
canonical programming. The last method was chosen because of its 
concise form. This method requires a minimum number of additions, 
subtractions, multiplications and storage elements. The time required 


for calculations is the limiting factor in the number of loops that 
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may be operated simultaneously. 


The D(z) transfer function shown in Figure 2.4 can be 


[3] 


expressed as the ratio of two polynomials in z 


E, (2) a +a 


he) 7 eae FE ae er ers (2.1) 


b toh ee te Og 
fe) 1 m 


D(z) = 


tA 


Where the coefficient of be represents the value of the weighting 
sequence at t = kT and k is an integer. The conditions for the D(z) 


transfer function to be physically realizable are that m > n and in 


dis) 


‘ =i 
the case where m = n, the denominator must not contain any factor z 


t 
and therefore b # 0. Equation (2.1) may be rewritten in the form 


0 
' Meo ae i 
a a, z a2 
SW + +e HO 
J t 
ae a to 
Dita) = ar se (2.2) 
by Zz bo Zz 
Ls : ee Me. a ; 
b b 
O fe) 


and, defining new coefficients, the relation becomes 


-1 =f 
E, (2) a, + az + se sak az 


Eid. cm in Tc Rita Gr ed Kea (2.3) 


-1 
1 i“ biz + 54 Me biz 


D(z) = 


A 


The condition for this equation to be physically realizable is that 
b # 0 unless a. # 0. Equation (2.3) is said to be of order n. 

Cross-multiplying equation (2.3) and taking the inverse z-transform 
yields 
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Figure 2.5 Realization of D(z) Transfer Function 
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' Solving for e, (t) yields 


y ae) (t - kT) - y bey C=) kT) (2.5) 
k=1 


Figure 2.5 shows the realization of the above equation. Since the 
output SG) depends upon present and past information of the input, 
it is advantageous to calculate the sum of the past information thus 
minimizing the number of calculations being performed when calculat- 
ing e,"(t). 


Referring to equation (2.5), we can calculate 
-) b,T (2.6) 


y (t) (2.7) 


M 
a 
a 


The above values are obtained before each sample period. 


At the sampling instant the calculation performed is 


e, (t) = ajfe, (t) - x (t)] +y¥ Ct) (2.8) 


The program accepts controller parameter definitions in 
the form given by equation (2.3). When a user requests a definition 
or redefinition, the program will request an order of controller in 
the range 1 <n < 5. The program will then request the necessary 


coefficients which are entered in engineering notation. 


2.5 PID Implementation 


The PID process controller is normally realized by analog 


’ gy Oy A oe 
pl eae eg ni 
: Shwe alla | 


ea | eas ae 
(25S) (T= 2)" 8yd_ <- (it = 3)" : 
faa oe 


% 


ais gomte no E8p9 | evods odg. Qo no. jesblbe: 
.sugat of3 to istacutguck i bas ; aRe 


evs notsemtoint te8q, 13, 20} mre 


| ae by Pel ‘ 7 “i 
’ 


~teluolso asdw bemrotzeq sated enotaetooles to odin addi 4 


Pat) 4 
cA 


| ' | 
sisinoiss aso sw ea s) notienpe oo ei 


ine 


(3 .S) 


(t.S) 


-bolisq eiqmse dos etnias bantetdo, ete ‘aulav sin 
ak bsiriotreq: notoaluotas sdt apse a 


dis ite cs = 


ak enokatatish TsJomoTE xelioxsaba. Siqsoon 
cbttiniteb »)sadnuins ony i “ ° mime 


18 


means. A digital computer may be programmed to generate the functions 
of differentiation and integration. The advantages of a digital 
computer controlling the process is that it can perform data logging 


and provide alarm indications with minor changes to the control loop. 


The digital PID controller is shown in Figure 2.6. 


Digital 
Computer 


Figure 2.6 PID Control System 


Proportional control is implemented by multiplying the 
x 
input ey (t) by the proportional gain constant S Integration of 


4 
the error signal ey (t) can be generated by 


* # * 
e,, (t) = e,, (t-T) “ Te, (t) (2.10) 


where en, (t) is the value of the integral, en, (t-1) is the value of 

the integral at the last sampling instant and T is the sampling rate. 
This method of control has the undesirable characteristic of being 
inaccurate for large step size, but this will not affect the integral 
control when the sampling period, T, is small compared to the major time 
constant of the controlled system. Since the process is closed-loop, 


the error will not increase without bound and in fact will tend towards 
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zero. 
Derivative control may be obtained by finding the slope 


of the line between the present value and the previous input sample 


or 
en 4 (t) =< [e,“(t) z e,” (t-T)] (2.11) 


x & 
where er (t-T) is the previous sample input and ey (t) is the present 
sample input. 


The output of the PID controller is given by 


* * * cs 
e, (t) = Se (t) + Keo, (t) + Keogh) Coie) 
where S = proportional gain constant 
K, = integral gain constant 
Ky = derivative gain constant 


The PID program requires T, ce K eT and K,/T to be specified for 
the control action desired. One output is computed and is checked 
for an overflow. If it is not in the range -10 V < ey (k) <.t10.V, 
the program will limit the value to the D/A converter's maximum out- 


put. Figure 2.7 shows a block diagram of the PID controller. 


Figure 2.7 Digital PID Controller 
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CHAPTER 3 


This chapter contains results obtained using all modes of 
the system program. All tests used the four A/D and D/A channels. 
The program can handle more processes but there are only four D/A 
channels and thus all results are shown for four systems running 
simultaneously. Recordings were obtained using a Brush 8-track pen 
recorder, with four tracks for the plant output, and four tracks for 


the control output. 


Sel D(z) Controllers With Common Sampling Rates 


A deadbeat controller was derived for the plant G(s) = 
1/s(st+1) from the state equations. Four plants were realized using 


standard analog techniques. The D(z) controller used was 


-1 
D(z) = Av nese , T = 1 second 


1+ .418z. 


The system response is shown in Figure 3.1 with the controller out- 
puts shown in Figure 3.2. 
Figures 3.3 and 3.4 indicate system response using a sampl- 


ing rate of 0.7 seconds for the following systems: 


r(t) = 1.5 Volts T = 0.7 seconds 
Figure System Controller 
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3.3-2 ey me aan ay) DL) Sea ee 
3.4-2 1+ 0.471z 


A\d wot yino sis aren tod aoe 
gntama ameteya su0i oe wipe sa 8 op: tie yA bas 
es a pe had > Se 

23g Soni4ee deuid s gates bombeado ar3W eo svleue mi ’ 


702% etosz3 xuot bas .3uqiuo scm sii 702) 


= (3)0 3osig eis toi baviaeb BBW astlorinog. 


gates bestisst stew see ES aw0% aii | 


21 


ides { ' \ : i 
eee Ge ee tees ty ei i| Pp iean 
ut 14 ea ‘ rt ‘ty | fect | Ma 
eet fee et ne re fa a EAE) ono Senne rae Mea 7 Rasa SIS. G, by Reed Nr ANE EOE 
[ON a | ; i a Te ti i H | Ain aah 
' coon fs , 4 feteeaesd a | ! 


Pe 


et (ee ir 


viper ii 
4 a = 
ey ! 


re ae ere ae tN eC ee - 
Pee Fir See 
S A i : t 
t : I Pas ‘ $ 
b i F ea 
, ' ‘ i 
+ + i 


ey 


fi 


Uden Tey ack Bla cttiet 
ocean eceey t+ watt tect 
' ati 1 


| Penis Zhuo aes ea tad 
bobo ne idee) ariesteir 


| ' 
' ' 
' | ' 

iif i! 


: Sees eo 
| 


are | ete orm eat 


I 


aaa ae rae aria 
! [ ; aie ' F ri ; i ' aa) eee mee | ' i i j ; 
Lf MOA SOPOUL. | | mya coro | ||| wmyancoroc | 


| 
UE ltssiak nee. i Pipe 8 kenemr at | | Lich Fier ete ane 1 Pri Sa! Been eee 


tte ft 


T= 1 sec. 


Response of systems with second order 


D(z) controllers. 


Figure 3.1 


‘ i se 
' 
es 
6 y 
-_ —s by g 
ate | J pe? ’ . 
ae on J: = ® 
™ > . 


on te pp ae he 
hoa f : 


Lt a de ge 
: } 


{ i 
pe hip apes 


i 
\ 

a 

; 

} 
a | 


; é 
f 
{ 


a 
\ 
ae @ 


ub ae | 
Per | 
ms At y ; 
‘a F —~. 
2) » * 
aii 
J 7 a a 
~4 é pore de 
h¢ OG 
J | 
* bape 


| 
} 
aa 
4 
| 
| 
| 


| 
v 

| 

| 


i 
} 


/ 


Hy 

. 

| 

P 
i 
| 

iA 


A eR, tm mG 


poe 


aad 


? ; : 
; ios ss 
ae t ar 
; i i L 
{ ; : 


i. 


ts = 7 
eh cert Beet Ee 
. | re. - 


4 i 
hy ee = . 2 ~~ 
a * 


( } 
] , ‘ 3 i 
ec i ee Rm alk 
H t 
e ali 
eh Se eo « 
és : mo : 


i i ¢ ; ‘ge 


; - i ps ms as 7 
Saha the SE BN, bg ee eae 
‘ 
. } 
‘ in 
—————< pe ). oo oe ae Ss 2 


{ \ 4 a8 j i ‘~~ 
ae? penne Se tm ey Mee, aed ee nt ar) 
: t : [ Le ta 


22 


ae ieee ay Seem) Sete 
owl igen 


SS eee 
fi t 


hs: ae ee SS a enn aceon fe 
ae ; 1 aan 4 
i { H 


atmo 


i 


Controller response for systems in Figure 3.1. 


a ed 


te as pref weer ad reac tartrate naam a 


i | ‘ 
ee 


Figure 3.2 


23 


CLS ia is CRLF eco 
Fife ee +i — sitter 3 Ps) 


T= 1 sec. 


Response of four differerit systems. 


Figure 3.3 
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Figure 3.4 
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Figure System Controller 
zi 
313.3 1 567 +1, Sez 
G(s) D(z) =-——————— 
4- : 
pas s(s + 0.8) 1+ 0.4542 2 
=I 
3 gf 1 35792000058 
i GAS) 39 ee aria ea aa bangs UL 
3.4-4 s(s + 2) Leute G82 


The above controllers were chosen to be deadbeat as this illustrates 
the most critical control functions. 

The next test involved controllers of various orders but 
with a constant sampling rate of 1.0 second. The plant used was the 


same for all cases. 


x 


r(t) = 1.5 Volts T = 1 second G(s) Sedat): 
r2eure Controller 
=a) 
3) 521 _ 1.58 - 0.582z 
4,621 D(z) = a) 
: 1 + 0.4182 
a =2 
3.5-2 .760 + .543z °° - .303z 
3.5-2 Dis a ed coe 
ies 1+ 0.721z ° + .217z 
3.5-3 815 - .3892-- + .88822 - .31472.> 
3 U2 ae =] -2 =3 
3.5- eee 17 ay Be e262 
a6 3) 6272 + Hehe i> $1280eL + eho  & 190s” 
D2) 5 = = =3 =4 
3.6-4 1+ .92° + .65z ~“ + .35z2 ~ + .095z 


The above controllers (except for the first case) were calculated with 


the requirement for a deadbeat response in a minimum number of sampl- 


ing periods removed. 
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Response of systems with 2,3 
order D(z) controllers. 


Figure 3.5 
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Controller response for systems in Figure 3.5. 


Figure 3.6 
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3.2 D(z) Controllers with Individual Sampling Rates 


The last test of D(z) controllers was performed using second- 
order controllers, but with different systems and different sampling 
rates. The systems are shown in Table 3.1. Again, deadbeat controllers 
were chosen as these provide the most severe test of the system. In 
actual practice, the input to the system will change gradually and the 
output will be less affected. 

A ramp input was not applied to the system since the analog 
computers did not have remote control of the OPERATE function. This 
is necessary in order that the start of the ramp he synchronized with 
the start of a sampling period. This feature was not placed in the 


operating system since it would not be used in the control laboratory. 


Ooo. LD. Controller 
A step input was applied to the PID control loop. The output 
of the three control actions is shown for a sampling rate of 1 second 


and 0.1 seconds in Figure 3.9 for a step input of 1 Volt. 


Figure Sampling Rate Control 
369-1 0.1 and 1.0 second A = 3.0 
359-2 0.1 and 1.0 second Ky = 0.1 
529-3 0.1 and 1.0 second Ka = 3.0 


Figure 3.10 shows the response of the proportional - 


Integral-Derivative outputs for the following closed loop systems: 
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Response of systems with sampling rates of 


0.2, 0.4,.0<7.,. and 2.0<seconds:. 


Figure 3.7 
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Figure 3.8 
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Response of PID control elements 


Figure 3.9 
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Response of PID control elements for the system 
T= 


1/ (2s+1) (10s+1). 


Figure 3.10 
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Figure 

3.10-1 Plant Output 

3210-2 Proportional Output 
3.10-3 Integral Output 
3.10-4 Derivative Output 


3.4 PID Controllers with Common Sampling Rates 


Using the same plant as above, and a sampling rate of 1.0 


second, the system response is shown for a 1.5 V step input. 
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The output of the controller in each case is the sum of all control 


actions. 


3.5 PID Controllers with Individual Sampling Rates 


Again, the same plant was used with a step input of 1.5 Volts. 
In this case the sampling rate was varied, indicating the effect of 
the sampling rate on the system response with control parameters fixed. 


Figure 3.12 shows the output response for the following cases: 
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Figure 3.13 shows the system outputs for the following 


Ks K, Ki 
3, OF45,; 6 
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The above results indicate that the PID controller can be 


useful where the constraints on the plant output are less rigid than 


for the case of the deadbeat controller. In addition, adjustment of 


the three gains can be varied until a satisfactory response is 


obtained. 


3.6 Calculator Mode 


sample calculations are presented in Figure 3.14. 


With the calculator routine running as a foreground program, 
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CHAPTER 4 


4.1 User Experience 


The controllers used in obtaining the test results in 
Chapter 3 were not affected by the delay (At) between obtaining a 
sample input, Cur and the output of a control, e,"(t). Table 4.2 
shows the delay, At, for the four system programs. This delay is a 
result of the calculations that must be performed in order to output 
a value which is dependent on the input. 

The time constant of the plant must be an order of magni- 
tude larger than the delay At. This is normally the case in physical 
processes where mechanical reaction time of a valve or pump is many 
orders of magnitude larger than the processor speed. 

When testing systems where a free integrator is present, 
the analog computer must be operational at all times with the input 
being applied via a switch. This allows zeroing of the integrator 
by its presence in the closed loop with no input r(t). 

Table 4.1 shows the maximum number of loops that can be 
run simultaneously for a sampling rate of 1.0 second. Using a sampl- 
ing rate of 0.1 seconds will necessarily reduce the number of controls 
loops for systems 0 and 2 to 1/10 the value shown in Table 4.1. 

Since all control parameters are entered into the execution 
table as the program is running, the analog computer should be in POT 
SET mode during the entry. The floating point routine will accept an 
indefinite string of characters but will only use the first seven of 


these. All inputs which are invalid will result in a ? being printed 
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and the user must either enter the value in correct form or request 
a new input. 

All mathematical calculations are checked for overflows, 
and in the event of an overflow, the value will be set to the 12-bit 
binary word representing +10 Volts or -10 Volts. Thus care must be 
taken in scaling all inputs to the system, if correct controller out- 


puts are expected. 


4.2 Conclusions 

The program presented in this thesis will allow the under- 
graduate engineer to become familiar with digital process control. 
The process can be considered to be of the continuous flow type as 
used in a chemical plant, a petroleum refinery, a cement plant or an 
electrical generation plant. Control of these processes refers to 
short-time control for providing stability and to long-time control 
for the purpose of maintaining some steady-state level of operation. 

Availability of this digital control program will enable 
the study of many varied and practical hybrid systems without program- 
ming in machine language. By modifying the floating point represent- 
ation of control parameters, it was possible to perform the calcula- 
tions for a control loop in the minimum possible time. Hence, a 


large amount of processor time remains and this will facilitate future 


program expansion. 
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APPENDIX A 


The following listing gives a brief description of the main 


subroutines used by the program in Appendix B. Flow charts for the 


complex subroutines follow the listing. 


IN 

OUT 
CRLF 
INTRP 
START 
RESTRT 


SAMPLE 


SYS1 


SYST1 


SERV 


CALDZ 


SYS 


COEF 
Pee 


FLPT 


MULT 


input a character from the teletype 

output a character to the teletype 

type a carriage return and line feed 

interrupt handler 

generates operating system 

restarts program without generating a new system 


saves all accumulators and carry before proceeding 
to control loops 


checks if time-out has occurred for entrance to 
D(z) control loops with common sampling rates 


calculates outputs for D(z) control loops with 
common sampling rates 


calculates outputs for D(z) control loops with 
individual sampling rates 


calculates values needed by the D(z) controllers 
for the next sample period 


dismisses sample interrupt and restores the pre- 
interrupt states of the accumulators and carry 


handles request for entry of D(z) parameters 
places D(z) parameters in appropriate storage area 


converts a floating point number to a modified form 
of floating point 


multiplies modified floating point number by an 
integer value 
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PIDO calculates PID control outputs for individual 
sampling rate loops 


PID) calculates PID control outputs for common sampling 
: rate loops 


TIME . converts a decimal integer to octal for clock timers 
STTL services a teletype interrupt 

MTYPE types a message stored in .TXT format 

OVFLO checks for overflows in control calculations 

END dismisses teletype interrupt service program 


CMODE a calculator routine using floating point arithmetic 
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Figure A.4  SYST1 - D(z) Control Loops with Common Sampling Rates 
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Figure A.7 PID1 - Common Sampling Rate PID Loop 
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Figure A.9 MULT - Multiply Routine 
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APPENDIX B 


sDIGITAL CONTROL SYSTEM 
3PKOGRAM IS SELF STARTING 


g00000 eLOG. | a 
GO000 BoRGR0 g 
00041 Zo9400 INTRP $ADDR OF INT 
BOOG2 O02046 STRT:  JMP @.STRT 3START PROGRAM 
90003 082047 RSTRT: JMP @eRSTRT 3RESTART PROGRAM 
3FOLLOWING LOCATIONS USED BY FLOATING POINT 
BOOOGA -LOC 4 s INTERPRETER 
200005 OCH gs s INITIALIZATION 
DLOOVE ~LOC 6 $TEMP WORD STORAGE 
BG0007 OO 4 7 $WRITEABLE AREA 
BO0B7 BO5300 WA: 5420-100 
3USED BY FLOATING POINT 
D2OB40 -LOC 48 
O2V40 900360 IN $INPUT A CHARACTER 
G2041 000366 OUT sOUTPUT A CHARACTER 
BBO0A4 -LOC 44 
@0044 G@1367 TTIIN: STTI sTTY SERVICE ROUTINE 
BQG45 BOGVOG RTCIN: @ 3ADDR OF SYSTEM SERVICE 
@0046 001056 .STRT: START sSTART ADDRESS 
@2047 G01307 .RSTRT: RESTRT 3RESTART ADDRESS 
B2O5@ OGG201 CMASK: 261 3CURRENT MASK 
G0051 00201 CMI: 201 sIDLE MASK 
@0052 GBG203 CMe: 203 3MASK OUT TTY 
00053 Q0G000 ORD1: B 3ASCII ORD OF CONTROLLER 
20054 GOGG0G ORD2: GB sOCTAL ORD OF CONTROLLER 
90055 177777 FLAG: -1 3TIME OUT FLAG 
99056 @0G407 -CLOC: MLOC-1 3ADDR OF D¢2Z) 
9Q@57 OGG497 «CL1: MLOC-1 3 PARAMETERS 
60668 080417 .CSTOR: MSTOR-1 3ADDR OF DCZ) DELAYS 
@2061 20417 -CSi: MSTOR-1 | 
@0062 177777 NDEV: ~-1 3# OF DEVICES ON LINE 
09063 177777 -CDEV: -1 3# OF DEVICES ON LINE 
@9064 002144 «MULT: MULT sMULTIPLY ROUTINE 
09065 900124 .TAB1: TABI-1 3TABLE ADDRESSES USED 
Q@0066 G20144 .TAB2: TAB2-1 3BY AUTOINCREMEMENT 
00067 GGG154 -TAB3: TAB3-1 3 REGISTERS 
00070 000164 .TAB4: TAB4-1 
@2071 920174 -TABS: TABS-1 
GOB7T2 GOGO CNUM: J 3DEVICE REQUESTING SERV 
90073 02240 -CMODE: CMODE-1 3CALC MODE 
@G074 999214 -MODE: MODE-1 
GO075 941616 -CEND: END sEND OF INTERRUPT 
@@076 GG1351 -MTYPE: MTYPE sMESSAGE TYPE ROUTINE 
90077 GG2065 -FLPT: FLPT sADDR OF FLOATING INPUT 
62180 GAGBEH -C6B: "BO 3ASCII CODES 
@G101 999061 -C6lz: “1 
G41G2 BOGG62 -C62: "2 
G0193 9OGG63 .C63: "3 
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TABS: 
TAB4: 
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MODEs 
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e BLK 
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3 CONSTANTS 
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3SUM BT'S TABLE 
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$TIMER TABLE 
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3 INPUT A CHARACTER FROM TTI 


CARRIAGE RETURN AND LINE FEED 


3CLOCK INTERRUPT? 
SERVICE IT 

TIARINTERRUPT 2? 

SYES> SOERVIGE “TTI 

3NO, UNIDENTIFIED 

3 INTERRUPT 

3LOAD DEVICE CODE INTO 


3YESs 
SNOs 


IN: NIOC TTI 
SKPDN Tt 
JMP e-1 
DIAC TTI 
LDA 1 »MASK 
AND 150 
3QUPUT A CHARACTER TO TTO 
OUT: SKPBZ TTO 
JMP o-1 
DOAS G»TTO 
JMP O53 
3TYPE A 
CRLF: STA 3s RTNG 
LDA OsCR 
JSR OUT 
LDA OsoLF 
JSR OUT 
JMP @RTN@ 
sINTERRUPT SERVICE ROUTINE 
«LOC AOD 
INTRP: SKPDZ RTC 
JMP @RTCIN 
SKPDZ TTI 
JMP @TTIIN 
JSR @oMTYPE 
MESS8& 
INTA @ 
HALT 


3ACC@)s5 AND HALT 


3 ADDRESS OF PARAMETER TABLES FOR 
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VARS+21 
VARS+42 
VARS+63 
VARS+1@4 
VARS+125 
VARS+146 
VARS+157 
VARS+226 
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OBGABS 
940513 
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082437 
824063 
G44062 
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26020 
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004430 
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O00773 
OOBA61 


G22023 
G2e29g2e4 
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3CONTROL LOOPS INCOMPLETE WHEN 
$IT IS TIME TO ENTER THEM AGAIN 


TMOUT: JSR @.MTYPE 
MESS9 S3INDICATE TIME-OUT 
HALT 3 AND HALT 
3COMMON SAMPLING RATE CLOCK TIMER 
Si ol’ NIOS RTC 3 RESET CLOCK 
DSZ DT 3 TIME FOR LOOP SERVICE? 
JMP e+5 3NOs RETURN 
STA ®»SAMP sYESs SAVE ACC®O) 
LDA OsN 3 RESET CLOCK COUNTER 
STA OsDT 
JMP SAMPLE+13 PROCEED TO SERVICE 
INTEN 3 ENABLE INTERRUPTS 
JMP eg 3 AND RETURN 
3 INDIVIDUAL SAMPLING RATE ENTERS HERE 
SYS@s NIOS RTC 3RESET CLOCK 
SAMPLE: STA @»SAMP 
ISZ FLAG 3LOOPS COMPLETED? 
JMP TMOUT 3NO»s ERROR! 
STA l»SAMP+13SAVE ACC 1) 
STA 2sSAMP+23 SAVE ACC 2) 
STA 32SAMP+33 SAVE ACC 3) 
MOVR 05@ 
STA Qs SAMP+43 SAVE CARRY 
LDA 0s@ 
STA @»SAMP+53 SAVE LOCATION @ 
LDA Qs CMASK 
STA OsSAMP+63 SAVE CURRENT MASK 
LDA Os CM2 
STA @sCMASK sMASK OUT TTY 
DOBS Os CPU 3AND TURN ON INTERRUPT 
3SYSTEM GENERATION WRITES INTO THIS LOCATION 
PUT 1: JMP @.13 
LDA l»eCDEV 3RESET NUMBER OF DEVICES 
STA l,sNDEV 30N LINE 
LDA 1»eTABS 3sGET TIMER ADDRESSES 
STA 1520 
LDA 15@20 
DSZ @2g 3LOOP TIME OUT? 
JSR IDLE! 3NO, IDLE THIS PERIOD 
JSR SERV 3YESs» SERVICE LOOP 
ISzZ NDEV 3FINISHED ALL LOOPS? 
JMP s=5 3NOs DO NEXT LOOP 
JMP UPDATE sYESs DO CALC FOR NEXT 
3’ SAMPLE PERIOD 
3 IDLE LOOP 
IDLEl: EDA O»@23 3 INCREMENT ALL 
LDA Qs@24 3POINTERS TO NEXT 
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3CONTROL LOOPS 


3INCR CHAN # 


3SET MODE=@, THEREFORE 
3NO CALC TO DO BEFORE 
3NEXT SAMPLE INPUT 
3TIMER TO EQUAL LENGTH 
30OF A SERVICE LOOP 


3 RETURN 

335 DECIMAL 
3RETURN ADDR 

3 ADDR OF SERV PROG 


3 INDIVIDUAL SAMPLING RATE D¢Z) CONTROLLER 


SERV: 


LDA 
DOAS 
DOB 
STA 
LDA 
STA 
LDA 
LDA 
LDA 
SKPDN 
JMP 
DICC 
SUB 
STA 
JSR 
LDA 
ADD 
JSK 
DOA 
LSz 
JMP 


lseE1 
1sADCV 
1sDACV 
39 -e1l2 
Qs@24 
Qs@24 
35@26 
1,@25 
232025 
ADCV 
e-l 
Q»sADCV 
350 
Bs@23 
@eMULT 
1,@27 
150 
OVFLO 
@»DACV 
°E1 
@.le 


3GET CHANNEL 
3START A/D 
3SET D/A CHAN 
3 SAVE RETURN 
SRE-INITIALICE 
3CLOCK TIMEOUT 
3GET SUM BT's 
3GET AG i 


3 A/D DONE? 
3NO, WAIT 
3YESs LOAD RESULT 


3 SAVE FOR UPDATE 
3AQ*CE1-SUM BT'S) 

3GET SUM AT'S — 

3 AT'S+AG*CE1-SUM BT'S) 
3CHECK FOR OVERFLOW 
3D/A CONTROL 

3INCR CHANNELS 

3 RETURN 


3STORAGE AREA FOR AC'S AND CARRY ON INTERRUPT 


SAMP: 


3ROUTINE UPDATES STORAGE 


- BLK 
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3ELEMENTS AND CALCULATES SUM BT'S AND SUM AT'S 


UPDATES: 


JSR 
LDA 
STA 
SZ 
PSZ 
DSZ 
JMP 


INIT 


1>-eCDEV 
1,NDEV 


« CLOC 


eCSTOR 


MODE 
0 t+2 


sINITIALIZE POINTERS 


3UPDATE REQUIRED? 
3NOs INCR POINTERS 
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TT 


09567 OBB405 JMP 05 3YES 
OQ570 B22G027 LDA Bs @27 
QO9571 @22026 LDA BW» O26 
Q@OG572 B22923 LDA B,@23 
@2573 888410 JMP +12 
@O574 626056 LDA 1,@.CLOC3GET ADDR OF PARAMETERS 
22575 O32060 LDA 2s@.CSTOR3GET ADDR OF DELAYS 
@8576 844030 STA 1532 
@8577 950031 STA 2531 
QOW6GB SWAAED JSR CALDZ 3CALC SUM AT'S, SUM BT'S 
6601 042927 STA O,@27 3 AND SAVE IN TABLES 3, 4 
QBCZ6G2 946926 STA 15@26 
@Q6G3 G18762 1SZ UPDATE+53INCR POINTER TO MODE 
08684 019062 1SZ NDEV 3 FINISHED? 
G@G6G5 OGH756 JMP UPDATE+33NO, DO NEXT LOOP 
Q@Z6G6 GH4436 JSR INIT 3YES 
@Q6G7 24115 LDA C2 3RESET CHAN # 
@8610 944447 STA i» OB) 
@0611 024432 LDA i,INSTR13 RESET INSTR IN 
@G612 944753 STA 1,UPDATE+53UPDATE +5 
8613 G24074 LDA 1,eMODE 
20614 220116 LDA Gs &0 | 
@8615 944920 STA 1520 3GET ADDR OF MODE TABLE 
@B616 30963 LDA 25eCDEV sPLACE 1°S IN TABLE 
@G617 G42620 STA B,@20 : 
0620 151404 INC 252,SZR 
QBW621 SGH776 JMP ~-2 
QG622 924061 LDA 1>eCS1 $3RESET POINTERS TO 
@0623 B44060 STA i»eCSTOR 3 PARAMETERS AND 
Q0624 G24057 LDA 1seCL1 sSTORAGE AREAS 
@0625 044056 STA 1seCLOC 
3ROUTINE TO DISMISS LOOP SERVICE INTERRUPT 
Q@0626 9608277 SYS: INTDS 3DISABLE INTERRUPT 
QG627 214855 DSZ FLAG 3RESET TIMEOUT FLAG 
QQV638 G20727 LDA Qs SAMP+6 
00631 PABO5O STA @»CMASK sREPLACE MASK 
Q@0632 062077 MSKO g 3MASK DEVICES 
8633 920722 LDA @»SAMP+4s REPLACE CARRY AND AC'S 
20634 1011080 MOVL 020 
98635 B2G714 LDA @» SAMP 
0636 @24714 LDA 1s SAMP+1 
@0637 038714 LDA 25 SAMP+2 
@0649 034714 LDA 3»SAMP+3 
Q@B641 @68177 INTEN 
@G642 002714 JMP @SAMP+5 
@0643 @14215 INSTR1i: DSZ MODE 
06644 924071 INIT: LDA l»eTABS sRESET AUTOINCREMENT 
QB645 B44G24 STA 1,24 3 REGISTERS 
QG646 B2407@ LDA 1,-TABA 
@Q647 B44023 STA i.23 
20659 G24065 LDA 1,-TAB! 
90651 G44925 STA 1,25 
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@5GG32 
024432 
G44430 
022031 
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3 ROUTINE CALCULATES SUM BT'S, 


STA 
LDA 
STA 
JMP 
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3A/Ds D/A CHAN # 


SUM AT'S AND 


3UPDATES STORAGE ELEMENTS 


CALDZ: 


eCC@s 
CCl: 


STA 
SUB 
STA 
STA 
LDA 
STA 
STA 
STA 
STA 
LDA 
DSZ 
JMP 
JMP 
LDA 
STA 
JMP 
LDA 
STA 
STA 
STA 
LDA 
STA 
LDA 
LDA 
LDA 
JSR 
LDA 
ADD 
JSR 
STA 
DSZ 
JMP 
LDA 
LDA 
LDA 
JSR 
LDA 
ADD 
JSR 
STA 
DSZ 
JMP 
LDA 
JMP 
i) 

1) 


3seCCG 
O5@ 
OseCC4 
BseCC5S 
1,@30 
1»-CCl 
1»eCC2 
2330 
2seCC3 
15,@31 
eCC1 
+2 
+4 
15,@31 
15833 
0-5 
1,@23 
15@33 
2531 
2532 
15-CC2 
15-CCl 
05631. 
25030 
15@30 
@eMULT 
15-CC4 
152 
OVFLO 
OseCC4 
eCCl 
e-il 
0,@32 
25030 
15@3@ 
@eMULT 
15-eCCS5S 
152 
OVFLO 
B»eCCS 
~CC2 
»-11 
I»+«CC4 
@.CCO 


3ZERO .-CC4 
3ZERO -CCS5S 
3GET ORD OF CONTROLLER 
3 AND SAVE 


3 SAVE ADDR OF STORAGE 
3 ELEMENTS 


3FINISHED SHIFTING? 
3NOs CONTINUE 

3YESs PROCEED 

3MOVE STORAGE ELEMENT 
3ONE DELAY LATER 

3 ANOTHER SHIFT? 

3PUT LATEST INPUT 

3IN FIRST DELAY ELEMENT 
3RESET 315 32 


3 RESET ORDER 


3$GET DELAY TN 
3GET BN 


3MULTIPLY 

3GET SUM 

3 ADDBN* TN 

3CHECK FOR OVERFLOW 
3 SAVE NEW SUM 
$FINISHED B'S*T'S? 
3NOs DO NEXT MULT 
3YESs GET DELAY TN 
3GET AN 


3MULTIPLY AN*TN 
3GET SUM 

3ADD AN*TN 

3CHECK FOR OVERFLOW 
3 SAVE NEW SUM 
SFINISHED A'S*T'S? 
3NO,s DO NEXT MULT 
3SUM BT'S IN ACC1) 
3 RETURN © 

3RETURN ADDR 

3ORDER OF CONTROLLER 


ESAT a2). AG ~~ 
Sak ate 
| | a) at eer nce ca 
® WARD ANG «GAs & : ta . chica 


Gia «* Ta Mud. eyT 


‘GPEA Mic s 


ATS 


ALUOJAD 


ALTUGH 


2TMIMZIS BaAGHNOTE CATA ‘ 
; 805446 AT2 ‘“:SGUA0 Seweeg. 
Gut U2. RAASET br 2800 
c KOG4 QHGS: ohOOsc8 At2 SEADOO Sage 
GUO. OF BSe 220148 ATe2 STARMAN OSs 508 
ABSICATVOOS. FO GRG Taar GEe «I! AQ’ genase haa 
. aVAzZ Ws IVa! ATe . BERRA eae 06 
SOQ» al ATe gebAap aa 
AROTe AO AGGA BVA. bus ATe EERHHAB rae 
CTVUAMS JS: 859648 ATS Top ges Bree " 
LER «! Ags) Le8ese 1Tee0° 
SOUITALRE Gena PMD & P39 SO boelG STaag 
WMETUOD OV ‘ay. aM Sheet. 2 ae 
CASA0RA .eayt 4ML AUAGHD AP Li 
TVAMSIS SVANOTS AVOMs 63.2 AGS reese “ere 20! 
AATAY YAMBU AMOy  ¢ Ate EeRee9 |araag 
STILE AaAHT OMAN Ge. tule EF Teae TTse a 
TUSMI TeNTad Tvs 62641 At ; —SS8SS0 aRYOR 
TUAMAUIS YAIAG TEAL wre EBB 0: AT2 EGO 1H 88 
SE 41S Teeads Pees ate 16ae20 saree 
| Se. |) — Seueee Bore 
AMGAO TAASH: SOD-vt AG Si ae Ge sovee 
| 1NGiee! AgTe BERD ABL epred 
“UT YA SAG T39 7% [68.6 AG ££0S86 aaras 9 


we 7 
vig a ee uy ‘ 


YIVIT A: = =TAUM«O AR PAVSOR If 
h MUG 2 Ae) 4 NOOval AGI 7 aShase sry 


1s MT eC 
WOLIRAVO. HOW .NODHD 
7 ‘MGE WAM SVAEL POD. 
7 “Se! THOR GQBHZ PUTA: 
THOM. FKav OG .OMe | theHs 
WT Yadad tao seg 
oe tee WA Twat @60nS8 


‘ oi 
; ak ' 


he fusaTsun 
Y a ae * ua ie T3I5% | 
a ie ae ogi \*) 
vi aa iV pee 


GUA t Ql 


7 4 bis 


; as 

y oy, 

: 7 | ‘? 
i a F 


i Ss 


} 


ACHSea TH 4 
REGd8Q OFF 


gages: pr ool 
2ESARG AITOO — 
(€SaQ00 @r7eo. 

TIREIG arrae | 
TOvOeO TETEO 
Séesse aspag 
WENSEN L 


4 
: a 
= 
pda 


oH 


tT? 


00736 
OB737 
QBTADO 
OO741 


06742 
00743 
OQVT4AA 
OBT4A5 
OB746 
OOTAT 
OB75G 
QO8751 
OB752 
QO753 
OO7TS54 
OB755 
OO756 
OBTS7T 
BWT60 
OO761 
00762 
04763 
OBT64 
OOBT65 
OB766 
OOTE7T 
OOTTO 
OO771 
OOTT2 
GO7TT3 
QO7TT4 
QB775S 
QOBQ7T76 
OOT7T7 
01000 
01001 
01002 
01083 
01004 
01905 
01006 
810067 
010190 
1911 
019012 
91013 
01014 
19815 
01016 
01017 
01020 
619021 


BOOOBBO 
UWOVBOD 
OBOGUO 
DOBBDO 


926123 
G40124 
924713 
M6512) 
66823 
036026 
026025 
032925 
063621 
OOOT77 
9626ea1 
162400 
BA2G23 
OBWEG64 
B2662T7 
123000 
Q™@4335 
061923 
018673 
910862 
Q607T56 
924115 
BA4667 
OB4422 
826434 
032435 
G44030 
050031 
PO4662 
BALV2T 
046026 
@1@425 
018426 
019062 
OOO766 
QGBAAGE 
924423 
44421 
G24417 
Q44415 
O@BGG614 
Q24963 
@4406e 
G24065 
B44Be5 
024066 
G4A4B26 
B24067 


eCC23 
“OCS? 
eCC4: 
eCC5s 


3 COMMON 


SYST 1,3 


INITe2: 
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3ORDER OF CONTROLLER 
3POINTER TO DELAYS 

2 ouM) BY Seth. S 

3SUM A*S«T'S 


SAMPLING RATE DC¢Z) CONTROLLER 


LDA 
STA 
LDA 
DOAS 
DOB 
LDA 
LDA 
LDA 
SKPDN 
JMP 
DICC 
SUB 
STA 
JSR 
LDA 
ADD 
JSR 
DOA 
1Sz 
Poe 
JMP 
LDA 
STA 
JSR 
LDA 
LDA 
STA 
STA 
JSR 
STA 
STA 
Sz 
| Sy a4 
1S2 
JMP 
JSR 
LDA 
STA 
LDA 
STA 
JMP 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 


OsN 
Q»DT 
1,-E1 
1, ADCV 
1,DACV 
33 @26 
1,@25 
232@25 
ADCV 
e-l 

05 ADCV 
35 
B5s@23 
@-MULT 
1,@27 
150 
OVFLO 
0»sDACV 
El 
NDEV 

SY OT bee 
15-C@ 
l>eEl 
INIT2 
1,@.T2 
25@.T2 
1530 
2231 
CALDZ 
Os@27 
15@26 
TO 

eT2 
NDEV 
e-12 
INIT2 
1,-T3 
loge le 
1,-Tl 
lseT@ 
SYo 

1, -CDEV 
1,NDEV 
1,-TABI 
1525 
1»e«TABe 
1526 
1s-TAB3 


3RE-INITIALIZE CLOCK 
3 COUNTER 

3GET CHAN ADDR 

3 START A/D 

3SET D/A CHAN 

3GET SUM BT°S 

3GET A@ ' 


3 A/D 
3NOs WAIT 

SYiS>s (LOAD RESULT &1 
3E1-SUM BT'S 

3SAVE FOR UPDATE 

3 A@*CE1-SUM BT'S) 

3GET SUM ATS ~ 

3 AT'S+AQ*CE1-SUM BT'S) 
3CHECK FOR OVERFLOW 
$3D/A CONTROL 

3 INCR CHANNELS 
SFUINLSHED, LOOPS? 

3NOs DO NEXT ONE 

SYES»s RE-INITIALIZE 

3 ALL POINTERS AND 

3 COUNTERS 


DONE? 


3CALC SUM AT'Ss SUM BT'S 
3AND SAVE IN TABLES 35 4 


SEENI SHED ALL LOOPS? 
3NO,s DO NEXT LOOP 
3YESs RESET POINTERS 


sDISMISS INTERRUPT 
3 RESET AUTOINCREMENT 
3 REGISTERS 
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01057 
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01066 
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01110 
01111 
O@1ile 
01113 
01114 
01115 
01116 
01117 


B4A4027 
B2LAB7TO 
944923 
QB1400 
000410 
OOB410 
OBGAL2G 
BOO42O 


ZO10B59 
OO2477 
177400 
060124 
177790 
02401 
OBG405 


068277 
BB6005 
268080 
QB6076 
082500 
G20772 
40774 
182400 
181400 
@24761 
44020 
230760 
Q42020 
151404 
Q00776 
B24755 
44020 
030754 
042020 
151404 
@B0776 
QG6076 
224070 
004360 
B24100 
030103 
142033 
186032 
04413 
122400 
24740 
123900 
49402 
804372 
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STA 
LDA 
STA 
JMP 
MLOC 
MLOC 
MSTOR 
MSTOR 
e EOT 


e LOC 
VAKS-1 
-400 
TAB1-1 
-18@ 
JMP 
JMP 


1,27 
1seTAB4A 
1523 
O23 


105@ 
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3 PARAMETER POINTER 


3DELAY ELEMENT POINTER 


3POINTER TO PARAMETER 

#- OF LOC TO LNI TI ALEEE 
3POINTER TO TABLES 

$1@@ LOCATIONS 


3ROUTINE TO GENERATE OPERATING SYSTEM 


START: 


INTDS 
FINI 
NIO 
JSR 
MESS7 
LDA 
STA 
SUB 
INC 
LDA 
STA 
LDA 
STA 
INC 
JMP 
LDA 
STA 
LDA 
STA 
INC 
JMP 
JSR 
MESS1 
JSK 
LDA 
LDA 
ADCZ# 
ADCZ# 
JSK 
SUB 
LDA 
ADD 
STA 
JSK 


4) 
@eMTYPE 


Os~eS6 
Os-451 
0s 

rn) 
lseSl 
1520 
2se02 
B»s@20 
esex oul 
e-2 

is 6S. 
1520 

ae ea4 
Gs@20 
23s2a95ER 
e-2 
@eMTYPE 


IN 

1s» -C6 
2seC63 
2202SNC 
O»s15SZC 
ERR1 
150 
PIeroD 
150 
Os251 
CRLF 


SINITIALIZE FLOATING PT. 
3MESSAGE INDICATING 

3 OPERATING PROCEDURES 
30F PROGRAM 

3TYPE MESSAGE ONLY ONCE 


sINITIALIZE PARAMETER 
3AND DATA STORAGE AREAS 


sINITIALIZE TABLES ON 
3PAGE @ 


3TYPE OF OPERATING 
SSYSTEM DESIRED? 


3 INVALID INPUT 
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01120 
21121 
01122 
01123 
91124 
01125 
01126 
01127 
01138 
01131 
01132 
01133 
01134 
21135 
01136 
01137 


01140 
01141 
01142 
01143 
01144 
01145 
01146 
01147 
01150 
01151 
01152 
01153 
01154 
1155 
01156 
01157 


01160 
01161 
01162 
01163 
01164 
01165 
01166 
@1167 
01170 
01171 
@11i7e 
Q1173 
01174 
091175 
01176 
01177 
91200 


B6009O 
081142 
081160 
001201 
091230 
1718900 
004356 
904372 
061370 
OBB4AA4 
002431 
000411 
004425 
QBO7TA2 
QBG644 
OQO707 


G2ed771 
G40045 
G2ea77T5 
B42573 
B2OSS4 
O42572 
O207T64 
842572 
920763 
942571 
020762 
042571 
086762 
820551 
042567 
OO@BATS 


094633 
B20T54 
042551 
920536 
B49045 
920535 
042551 
020743 
G942551 
920527 
B4255G 
0942551 
G42551 
942551 
920525 
042544 
0BB454 


NIO a) 3JMP @ TO GENERATE 
GEN@ sDESIRED SYSTEM 
GEN1 
GEN2 
GEN3 
ERR1: MOV 352 sINVALID INPUT 
JSR ERR+1 sTYPE ? AND ASK 
JSR CRLF 3 AGAIN 
JMP -1052 
GOO: SYS@ 
-GO1: 002431 3JMP @.FLI1 
°GO32 JMP +11 
-GO43 Q04425 3 JSR TIME 
-GOS: SYST 1 
~GO63 INIT 
~GO7? OQO7B7 3 JMP PIT6+1 
3GENERATE INDIVIDUAL SAMPLING RATE DC2) 
3 CONTROLLERS 
GENG: LDA Ds GOO 
STA @,sRTCIN 
LDA Ose GO7 
STA @»@PUT 
LDA Qs» NOP 
STA @s@PUT+1 
LDA Os -G@1 
STA O,@PUT+3 
LDA Os» -GO3 
STA @,@PUT+4 
LDA O»-GO4 
STA @»@PUT+6 
JSR @.GH6 
LDA Os G22 
STA Q»s@PUT+7 
JMP GEND 
3GENERATE COMMON SAMPLING RATE D(Z) CONTROLLERS 
GEN13 JSR INIT2 
LDA O»s~-GO5 
STA 0»@-G31 
LDA Os-G12 
STA B»RTCIN 
LDA OseGll 
STA @,@PUT+1 
LDA B»-GQl 
STA @»s@PUT+3 
LDA @»~eNOP 
STA @,@PUT+4 
STA Q»@PUT+6 
STA Qs @PUT+7 
STA Q»s@PUT+10 
LDA O»s-G14 
STA @»@PUT+5 
JMP GEND 
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id a 


3GENEKATE INDIVIDUAL SAMPLING RATE PID 
3 CONTROLLERS 


Y1i2eV01 O20738 GEN2: LDA Os- Gd 
YVl2V2 B49045 STA WstkTCIN 
01203 020521 LDA O0»5-G2d 
01204 042527 STA 0s@-G31 
01285 020513 LDA 6s -NOP 
01206 842530 STA 0, @PUT 
91207 842534 STA O@s@PUT+5 
01210 920512 LDA O»-Gil 
01211 842526 STA O»s@PUT+1 
Glelte @20513 LDA O,5-Gel 
91213 842525 STA @,@PUT+2 
01214 020514 LDA QseGe4 
01215 842524 STA O,s@PUT+3 
91216 020510 LDA O»s-G2e 
01217 042523 STA Os@PUT+4 
01220 G@20714 LDA B»s-GO4 
1221 842523 STA Os@PUT+6 
01222 G20504 LDA O»-G22 
01223 942522 STA Os @PUT+7 
G1224 8920510 LDA O5-G32 
01225 042521 STA 0»@PUT+19 
01226 086503 JSR @eG2s5 
01227 BBO425 JMP GEND 
3GENERATE COMMON SAMPLING RATE PID CONTROLLERS 
91230 820471 GENS: LDA Os-eG10 
01231 G40045 STA O»RTCIN 
01232 B82058G LDA Bs ~-G3G 
01233 BZA2500 STA Os@-G3l 
91234 820466 LDA Os-Gll 
01235 642502 STA O,@PUT+1 
01236 O204T77 LDA O0s+G33 
01237 842501 STA O,@PUT+2 
G1240 B2047B LDA O» G24 
01241 942580 STA @,@PUT+3 
01242 920456 LDA QBs ~eNOP 
01243 842473 STA @,»,@PUT 
91244 942476 STA O,@PUT+4 
91245 G42476 STA O,@PUT+5S 
@1246 B424T77 STA @s@PUT+7 
Q1247 B2G460 LDA O»seGe3 
91250 042474 STA @s@PUT+6 
91251 8020463 LDA B»s-eG3e 
G@1252 GA42474 STA O,@PUT+19@ 
81253 806456 JSK @.Ge5s 
3SYSTEM GENERATEDs, COMPLETE INFORMATION 
91254 @894475 GEND: JSR MTY PE 3REQUEST NUMBER OF 
1255 BG404O MESS@ sDEVICES ON-LINE 
01256 884369 JSR IN 
01257 824191 LDA 1» -C6l1 


@1260 030104 LDA 2s -C64 
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01631 
91632 
01633 
01634 
91635 
91636 
01637 
01640 
01641 
01642 
81643 
01644 
1645 
01646 
01647 
01658 


042021 
GA6621 
GB4372 
GBQ6B76 
OVA2B7 
OB6077 
942021 
946021 
102400 
942021 
042021 
094372 
QB6076 
084213 
OB6077 
942021 
QG46021 
QBWB4B5 
DOBVOD 
COOO71 
G@BBB12 
QBOBT7 


OBA372 
G602T7 
020613 
B40050 
@62077 
060210 
Q60177 
082685 


20100 
G40503 
034056 
175400 
G24072 
137000 
G24054 
031400 
132400 
132400 
132400 
132406 
158008 
151400 
151400 
150000 
ASOG22 
030065 
G24072 


eAls 

eC71s 
-C12s3 
QUES: 


STA 
STA 
JSR 
JSR 
MESS6+3 
JSR 
STA 
STA 
SUB 
STA 
STA 
JSR 
JSR 
MESS6+7 
JSR 
STA 
STA 
JMP 


Li2 
77 


Os@21 
15@2l 
CRLF 


@eMTYPE 


@eFLPT 
O,5@21 
15@21 
ee) 
Qs@21 
Os@21 
CRLF 


@eMTYPE 


@eFLPT 
Os@2l 
1,@21 
END 
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3 SAVE CHARACTERISTIC 
3SAVE MANTISSA 


SDEFINE KI 
3SAVE CHARACTERISTIC 
3SAVE MANTISSA 


ZERO) Bet CK=1) 
$ZERO El 


3DEFINE KD 

$SAVE CHARACTERISTIC 
3 SAVE MANTISSA 

3END OF INTERRUPT 
sOCTAL TIME SUM 


3RESTORE PRE-INTERRUPT STATE OF 
3 PROGRAM AFTER TTI INTERRUPT 


END 


JSR 
INTDS 
LDA 
STA 
MSKO 
NIOC 
INTEN 
JMP 


3D¢Z) PARAMETER 


COEF: 


LDA 
STA 
LDA 
INC 
LDA 
ADD 
LDA 
LDA 
SUB 
SUB 
SUB 
SUB 
COM 
INC 
INC 
COM 
SPA 
LDA 
LDA 


CRLF 


O»sTISAV+1 


@»CMASK 


1) 
TTI 


@TISAV 


DEFINITION AND STORAGE 


Os -eC60 
OseK4 


3»-«CLOC 


353 
1»sCNUM 
153 
1,0RD2 
23853 
1,2 
152 
1,2 
152 
222 
222 
222 
222 
2222 


2seTABI 


15CNUM 


3POINTER TO PARAMETER 
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ttt 
01651 
01652 
091653 
01654 
01655 
01656 
@1657 
01660 
01661 
01662 
01663 
91664 
@1665 
01666 
01667 
1670 
01671 
01672 
01673 
Q1674 
91675 
01676 
O1677 
01700 
1701 
01702 
01703 
01704 
01705 
01706 
Q1707 
21710 
Q1711 
Q\71le 
01713 
01714 
01715 
01716 
O1717 
01720 
01721 
1722 
01723 
01724 
Q1725 
01726 
91727 
91736 
01731 
91732 
91733 
01734 
91735 
01736 


133000 
133000 
050021 
Q2B457 
084366 
20457 
004366 
020452 
004366 
020454 
084366 
020122 
004366 
010444 
QQ60T7 
@190437 
BOO4O4 
842021 
946021 
OBB403 
Q42G022 
046022 
GB4372 
020053 
924431 
122434 
080751 
920425 
G24426 
106414 
GBO4G6 
920424 
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020101 
040417 
OOBTAD 
Q24054 
046022 
Geg4Ale 
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GBGETO 
LTT. 
177777 
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BOOGEG 
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PPT 1s 


eK1: 
Kes 
eK33 
eK43 
eKS8 
eK68 
eK7$ 
eKSs 


ADD 


JSR 
LDA 
JSR 
LDA 


1Sz 


SUBZ# 


SUB# 


JMP 
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101 
"oO 
TAN 
"RB 

be  ¢ 

is 2) 
-EOT 


1,2 
1s2 
2ee2i 
OseKS 
OUT 
OseKT 
OUT 
Os eK4 
OUT 
O»eK8 
OUT 
B»sEQ 
OUT 
°K4 
@eFLPT 
eK1 
e+4 
O,@21 
1,@21 
e+3 
O»s@22 
1s@e2e2 
CRLF 
@»sORD1 
ls eK4 


1»s@sSZR 
COEF+26 


OseK3 
lseKS 


Qx15SZR 


e+6 
O»seK6 
OseKS 
O»eC6l 
O»seK4 


COEF+26 


1,O0RD2 
1,@22 
Ox eK3 
OsxeKS 
Os ~eC60 
QseK4 
l»eKe 
ls eKl 
CRLF 
END 
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3POINTER TO A@ 
3 STORAGE AREA 


3 INPUT A NUMBER 


3SAVE AQ CHARACTERISTIC 
3SAVE AS MANTISSA 


3 SAVE PARM CHARACTERISTIC 
3SAVE PARM MANTISSA 


3A"S COMPLETE? 
3NOs GET NEXT A 
3YES 


3HAVE B'S BEEN ENTERED 
3YESs EXIT ROUTINE 
3NO»s GET THEM 


3SAVE ORDER OF 
3CONTROLLER IN PARM 
3 STORAGE AREA 

3 RESET CONSTANTS 


3END OF TTI INTERRUPT 
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001750 e LOC 1750 
3PID IDLE LOOP FOR INDIVIDUAL 
a¢ 3 SAMPLING RATES 
@175@ 024501 IDLE@: LDA 1seP6 3INCR REGISTERS TO POINT 


91751 Oe22922 LDA O0s@22 3TO NEXT PID LOOP 
01752 125404 INC ls1isSZR 
91753 G@B0776 JMP 07-2 
@1754 G22824 LDA Os@24 
Q@1755 822024 LDA 0»@24 
@1756 019466 ISZ «Pil S3INCR A/Ds D/A CHAN # 
@1757 824473 LDA lseP7 3WAIT EQUAL TIME OF 
1768 125404 INC lsisSZR sCONTROL CALCULATION 
91761 @O0777 JMP emit 
01762 O@BO456 JMP PUT2-1 $3$RETURN 

3 INDIVIDUAL SAMMPLING RATE PID CONTROL LOOPS 
@1763 924063 PID@: LDA 1seCDEV 3RESET # OF DEVICES 
01764 844062 STA 1 »NDEV 
91765 624071 LDA l»-eTABS 
21766 844020 STA 1520 
91767 B26G82o LDA 1,@20 
91770 01682e DSZ e280 3LOOP TIME-OUT? 
01771 @BATS7T JSR IDLE@ 3NO» IDEE THIS. PASS 
G1772 B26024 LDA 12@24 3YESs RESET CLOCK 
01773 G@46024 STA 15@24 


3PID CALCULATOR USED BY COMMON AND INDIVIDUAL 
3SAMPLING RATE LOOPS 


01774 @2445@ PIDI1: LDA lsePl 3GET CHAN # 

1775 865121 DOAS lsADCV sSTART A/D 

01776 966023 DOB 1»sDACV $SET D/A CHAN # 

91777 @26022 LDA 1,@22 3GET KP 

G2009 @3282e2 LDA 22@22 

G2091 863621 SKPDN ADCV 3A/D DONE? 

G2d02 BOOTT7T JMP ofl 3NO»s WAIT 

2003 962621 DICC @sADCV s$YES»s LOAD RESULT, E1 

B2004 B40447 STA OsE1 3SAVE INPUT 

02085 006064 JSR @eMULT sP2=KP*E1 

G2006 B40437 STA QseP2 3SAVE RESULT 

G2007T B2G444 LDA OsE1 

02010 G26022 LDA 1,@22 3GET KI*T 

92011 B832022 LDA 22@22 

92012 GB6064 JSR @eMULT sKI*T*E1 

2013 G2602e2 LDA 1,@22 3GET KI*EeICK-1) 
LDA 3222 

#8015 qa0n ADD 15@ 3 P3=KI¥*¥T*EL+KI*E1CK-1) 

2016 841400 STA O23 3SAVE SUM 

G2617 B40427 STA DseP3 3SAVE SUM 

2920 G2602e2 LDA 1,@22 3GEL E1CK=1) 

G2021 B2043e LDA OsE1 3GET E1 

92022 941401 STA Bs1l53 3E10K-1)=E1 

2023 122400 SUB 150 3D=E1(K-1)-E1 

G2024 826022 LDA 1,@22 3GET KD/T 


G2025 932022 LDA 2s@22 


_ BX 


JaudInia a S401. ate | 


i, 
OF anatepoan: AEA 
W003 GPS TAAIUN GTe 


TATOS 


ViAHD ANG GNA 
26 aMiT JAUGR 
4BOITAIUIIAD 


‘ AROMLA 


TTAWt 


AOA THEO & 


wart ‘aewl! 
nL 
SaTAR PUIUaMASL ; vie | 
AGU 299N07 Dense BereD 
Agu | “SSeSsn 1erre 


ow al 
SS0 <8 


WESulat OME - pOSeST Sete 
Cae Shi trees cette 

A898 Ags ASOSSE KATH 

. 8628 AG’ ASOsss ce ert 
‘PSs xe] 800018 derTe 

Thies AO EvAASO Terk 
NES hail OWd, sonens ee Wi:) 

b~4 “ML i 


Trage rer 


VARUTAA f~SFUS WL a@aoge $a) 

. sy an 

2900.) JOATVOO Gi9 3TAR PALISMMAL JAUGIV paws woe 
CaQIVSG AD & Ta2ant VaGowet AGS s8019% SBOASB : cavite 
VGGh ct Ate , Sa0saG AaTIO 

[HATS «tT AGA pes gaye 

Sek | A¥e SGAay aorre 

WSHel AGL osasso po eted 

T TU MET Qood: Hse SOG aseaie Qrrig 

2@A9 2IHT BiG) .Ouy | OFaaT He: Tevaeg ITTIe 
HIGJ aan «day LSBet Adu  ASHBSR Sx FE 
AS@.! Ags) 6S$0aso EFT TO 

‘ae ¥ rg ; , ee 

JAUGIVIGUI, QUA BOMMOO Ya seis HOTALUSZAD Ids i a 


* WAKO. Tae 
CNA 
* WAKO ANG Tae 
Se anol ty 


SayOd. d\Ae 
TIAW ,OUE. 
Ia «TIUCAR UAW eG a Y< 


TUSHI BVAL’ 


LAtANaASS ¢ 
TIER AVA 


eID 


TAD: 


be TLE. 
Ci-WNISS#IN Tabs 


Cie woES4 L+ 1 Rae Ts Diee4 \ 


. | MLgdVAz 
a ae Mue. aUADt 
Ch-H) IG TBar 


as 


[S-¢f-N) ng=ae 
TANG. T 


THAT 
c) MOAGE? 


“tS Tae. . 
[A=¢l-mLar 


we te 


@400.! ATAA DALSSMAZt | ‘J 
a: ae Ag’ -31GI9 Re : 


VadAal epoca | 
) Oe. 
we BT (on GROET) * 
SSags if Aan ite 
‘ VGA. nate - ysacag. 
i fms, ME | 
VOdAL® Sart 
f3<8 are 
Foes Sr Hel 
Qs” Ate” 
Mee apt add 


SSO t = 


Ste TF. 
Bees pes ver 


“sagas Bi y 


:  TPAg ED se08e 
ee Q 

ha 
istic feaae 
- SSes0- Hoss 
Bite 


Th 


Sseayg 
TAM et 
' SSeet 


rer 


TTT 


02026 
WG2edeT 
02030 
02031 
G92032 
02033 
92034 
02035 
02036 
G2037 
@2040 
02041 
92042 
B2043 
02044 
92B45 
G2046 
OG2047 
92050 
G2051 
B2gse 
92053 


OG2054 
2055 
92056 
QG2057 
J2060 
02061 
QW2062 
02063 
02064 


G2065 
02066 
Q2067T 
QZ20TO 
92071 
G20T2 
92073 
B2eOT4 
G2ed7s 
B2076 
02077 
92100 
92191 
92102 
92183 
G21B4 
92105 
021086 
G2197 


OG6064 
024416 
123000 
984335 
101020 
924413 
123080 
084335 
961123 
Q@10405 
G10062 
098733 
QBAALe2 
QB24B5 
QOOBDO 
BOBWOD 
BBOOWGOD 
GG2737 
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17777@ 
IG I7TA 
QOOOBO 


820773 
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1G240@ 
946763 
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940024 
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O54446 
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1290009 
BABAS2 
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O300G07 
V219G2 
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186405 
OBBAGE 
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028110 
904366 
BDA437T2 
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3KD*D 
3GET P2 

3 KP*E1+KD*D 

3CHECK FOR OVERFLOW 
3ZERO CARRY 
3CONTROL OUTPUT, E2 
3CHECK FOR OVERFLOW 

3D/A CONTROL 

3INCR CHAN # 

3ALL LOOPS COMPLETE? 

3NO» DO NEXT LOOP 

3YES, RESET POINTERS 
3EXIT CONTROL INTERRUPT 
3A/Ds D/A CHANNEL # 
3KP*E1 

$KI*XT*E1+KI *E21 

3ADDR OF PARMETERS 

3EXIT INTERRUPT 

3# OF AUTOINCREMENTS/LOOP 
3TIME FOR PID CONTROL 

3 INPUT £1 


3MODIFICATION OF FLOATING POINT INPUT 


FLPT: 


PUT33 


3SAVE RETURN 

sENTERs GET FLOATING 
sPOINT INPUT AND SAVE 
332=-BIT_ RESULT 

3AND EXIT 


3GET BREAK CHARACTER 


S15 °1T A CR? 
3YES 

3NO3 TYRE? 
3AND REQUEST 
3A NEW INPUT 


3 INPUT O-Ke 
3GET 16-BITS 
3SAVE SIGN 


JSR @eMULT 
LDA 1»5eP2 
ADD 152 
JSR OVFLO 
MOVZ OsB 
LDA 1seP3 
ADD 150 
JSR OVFLO 
DOAS QsDACV 
LiSsz eP1 
lisZ NDEV 
JMP PID1 
WSR INITS 
JMP @.PS 

i) 

i) 

i) 

PARM-1 

SYS 

-1@ 

-1 

i) 

POINTERS AND CONSTANTS 
LDA BseP4 
STA 0,22 
LDA Os eCDEV 
STA Bs»NDEV 
SUB OsO 
STA Os ePl 
LDA O»seTABS 
STA 0224 
JMP G53 
STA 3sRTN2 
FETR 

FDFC 4) 

FSTA BsNUM 
FEXT 

LDA 2eWA 
LDA Bs2s02 
LDA 15CR 
SUB B»sisSNR 
JMP +6 
DSZ eK4 
LDA Bs eCTT7 
JSR OUT 
JSR CRLF 
JMP @.FL1 
LDA 2sMSK2 
LDA OsNUM 
MOVZL Gs@ 
AND 
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OBTBOD 
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DOGABE 
125113 
100001 
100001 
BBWAVe 
GBW4O3 
G34452 
954452 
962001 
973301 
VASPS5 
GBB436 
030451 
G60401 
133404 


ERR2s 


RTN23 
MSKe2: 
MSK33: 
MSK43 
MSKS53 
MSK63 
NUM3 


JMP 
LDA 
ANDS 
JMP 
LDA 
ANDS 
LDA 
ADD 
MOVR 
LDA 
LDA 
LDA 
ANDS 
COM 
ANDS 
ADD 
JMP 
LDA 
JMP 
7) 


COEF+26 


B40B00 
OBTIBBA 
077080 
177600 
177402 
e BLK 


4 
2eaMSK3 
250 
et+5 
2sMSK4 
220 
2e2MSK5 
220 
Ds@ 
2esMSK6 
1 »NUM 
3sNUM+1 
223 
222 
21 
3a1 
@RTN2 
3seRTN2 
-3593 
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3NO 
3YESs PLACE POWER 
3RIGHTMOST BITS 


3SAVE POWER IN 
3RIGHTMOST BITS AND 
SRL ELOUN. dS TO SLERT 


3REPLACE SIGN 


3GET 32-BIT NUMBER AND 
3 FORM 16-BIT MANTISSA 
3RETURN WITH CHARACTER- 
3ISTIC IN ACC8)s AND 
s3MANTISSA IN ACC1) 


3REQe NEW INPUT 
3RETURN ADDR 


3REQ@-e NEW DCZ) INPUT 
3MASKS FOR MODIFICATION 


332-BIT FLOATING PTe # 


3ROUTINE MULTIPLIES CHARACTERISTIC IN ACCQO)s5 


s3MANTISSA IN ACC2) 
3AND RESULT IS RETURNED IN ACC@)- 
AND IF IT OCCURS THE 


31S MADE FOR OVERFLOW, 


3RESULT IS REPLACED WITH MAX. 


MULT $ 


CMULT : 


MULT 1: 


STA 
NIOC 
MO VL# 
JMP 
MOVL# 
JMP 
JMP 
MOVL# 
COM 
COM 
JMP 
JMP 
LDA 
STA 
DOB 
DOCP 
MOVZL# 
JMP 
LDA 
DIA 
AND 


35RTN6 
MDV 


QsBxS2C 


+4 


1515SNC 


MULT 1 
CMULT 


1s15SNC 
BsOsSKP 
@sO@sSKP 


CMULT 
MULT 1 
35CR1 
35CR2 
OsMDV 
2sMDV 


1515S5NR 


FIN-1 
25MSK 1 
@»MDV 


1s2sSZR 


BY 16-BIT NUMBER IN ACC1) 


A CHECK 


D/A VALUE 
3SAVE RETURN 


3ACC@) +VE? 

3NO 

3YESs, ACC1) +VE? 

SY ESs “MULTIPLY 

3NOs MULTs COMPL- RESULT 
3ACC1) +VE? 
3YES» COMPL- 
3NOs MULTIPLY 


RESULT 


3SET COMPLIMENT 

3 FLAG 

3LOAD MULTIPLIER 

3AND START 

sANY SHIFTS REQ*D? 

3NOs GET RESULT 

SYESs GET FIRST 16-BiTs 
30F MULTIPLICATION 
3SHIFT RIGHT? 
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92171 
@2172 
G2173 
Q2174 
92175 
02176 
02177 
B22oo 
02201 
G2e20e2 
G2203 
B2204 
G2205 
G2206 
G2egT7 
92216 
Geeii 
G2elea 
92213 
92214 
Cer ld 
92216 
92217 
92226 
G2221 
Vecee 
B2223 
G2224 
Ge2ee5 
92226 
Beeel 
92230 
92231 
Ge2e3e 
02233 
G2234 
92235 
92236 
G2237 
B224G 


2241 
gee4e 
G2243 
G2244 
G2245 
Q2246 
D22A4T7 
G2256 
geesi 
Geese 
G2253 


009421 
O54444 
125180 
125226 
124499 
071401 
934435 
1511090 
101100 
175494 
00775 
934431 
163414 
169000 
125494 
QOO767 
009413 
1251990 
125249 
191220 
1912290 
191220 
101220 
125404 
090773 
099402 
060481 
191920 
319406 
002403 
1090020 
002401 
BOGHAO 
La7797 
QBBBO 
177774 
174080 
02000 
G00020 
096005 


QDA3Te 
OBA372 
B20475 
0B43 66 
820515 
GGA366 
QBEGGA 
129000 
109090 
B38G07 
G21002 
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RTN6? 
CRit 
CR2: 
K1s 
OFLOs 
RTN13 
MSK13s 


3CALCULATOR MODEs 


JMP 
STA 
MOVL 
MOVZR 
NEG 
DIB 
LDA 
MOVL 
MO VL 
INC 
JMP 
LDA 
AND # 
COM 
INC 
JMP 
JMP 
MO VL 
MOVOR 
MOVZR 
MOVZR 
MOVZR 
MOVZR 
INC 
JMP 
JMP 
DIA 
MOVZ 
ISZ 
JMP 
comz 
JMP 
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-1 
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-4 
174000 
Q 
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FINI 


RT 
32RTN1 
1,1 
151 
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2sMDV 
339K 1 
2392 
a) 


3935SZR 
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1515SZR 


e-11 
FIN 
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FIN 
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Os@ 
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3YES 


3NOs SHIFT LEFT 


3GET LAST 16-BITS 


sSHIFTED 4 TIMES? 
3NO 
3YESs OVERFLOW? 
REPLACE WITH MAX. 
ANY MORE SHIFTS 
DO THEM 

FINISHED 


SYESs 
3sNO>5 
3YES>s 
3NO> 


3SHIFT RIGHT 4 TIMES 


3 FINISHED? 
3NO 
3YES 


3ZERO CARRY 

3 COMPLIMENT RESULT? 
3NOs RETURN 
3YESs COMPLe 
3 RETURN 


RESULT 


3COMLIMENT FLAG 


34 SHIFTS 
3OVERFLOW CHECK 


PERFORMS +5s-5*sAND / 


3USING FLOATING POINT INSTRUCTIONS 


CMODEs: 


JSR 
JSR 
LDA 
JSR 
LDA 
JSR 
FETR 
FDFC 
FEXT 
LDA 
LDA 


CRLF 
CRLF 
Bae FFA 
OUT 
@»sSP 
OUT 
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2sWA 
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3 CALCULATOR MODE 


3GET A NUMBER 


3GET BREAK CHARACTER 
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G2323 
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02336 
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92341 
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106415 
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QGBEGGA4 
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BBBTOE 
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OBOBST 
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e ADD: 


e SUB: 


eMPYs: 


oF FOS 
eFF 18 
oF F232 
eFFS33 
oFFA: 
eo FF 63 


LDA 
SUB# 
JMP 


SUB# 


SUB# 
JMP 
LDA 
SUB# 


SUB# 
JMP 


1seFFQG 
Qs15sSNR 
eMPY 
lse FFI 
Os1isSNR 
«DIV 
lseFF2 
Os15SNR 
e ADD 
lyse FFS 
OsilsSNR 
e SUB 
lseFF6 
Ox1sSZR 
eot+1]1 
OsoFF4 
OUT 
OseFF4 
OUT 
CRLF 
CRLF 
CRLF 
@.CEND 
CRLF 
OseC7T7 
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3MULTIPLICATION? 


3DIVISION? 


3 ADDITION? 


3 SUBTRACTION? 


3EXIT CALCULATOR? 
3 INPUT ERROR 
SINDICATE EXIT OF 
3CALCULATOR MODE 


3 AND RETURN 
s3TYPE ? AND WAIT 
3FOR ANOTHER INPUT 


3 ADD 


3 SUBTRACT 


yMUTETT PLY 


3DIVIDE 


3TYPE CRLFs = 


3TYPE RESULT 


3DO NEXT CALCULATION 
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3TYPE A CARRIAGE RETURN AND EQUAL SIGN 


eF3s 


SP: 


@ 


STA 
LDA 
LDA 
LDA 
SUB# 
JSR 
LDA 
JSR 
LDA 
JSR 
LDA 
JSR 
LDA 
JSR 
JMP 


Sse FFT 
2sWA 
Bs2s2 
1,CR 
@s15SNR 
CRLF 
BsSP 
OUT 
@»SP 
OUT 
O»sEQ 
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3PARAMETER STORAGE AREA 


VARS $ 
STOR: 
PARM: 


e LOC 


e BLK 


e BLK 


e BLK 
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3DCZ) CONTROLLERS 
3D€CZ) DELAY STORAGE 
3PID CONTROLLERS 


3MESSAGES ARE STORED IN LOCATIONS 2580 - 4405 
3FOLLOWING IS A BRIEF DESCRIPTION OF THEM 


3AND THE LOCATIONS THEY OCCUPY. 


IS 


sDESTROYED WHEN OPERATING SYSTEM IS FORMED 


3MESS7 
3MESS@ 
3MESS1 
3MESS2 
3MESS3 
3MESS4 
3MESSS 
3MESS6 
3MESS8 
3MESS9 
3MESSA 


DESCRIBES OPERATING SYSTEM 
# OF DEVICES ON-LINE 

# OF DESIRED SYSTEM 

SYSTEM INITIALIZED 


COMPUTER # 


ORDER OF D¢€Z) CONTROLLER 


SAMPLING RATE 
PID PARAMETERS 


INVALID INTERRUPT 


TIME-OUT 
RTC FREQUENCY 
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*** DIGITAL CONTROL SYSTEM xxx 


SYSTEM PROGRAMS 

@ - INDIVIDUAL SAMPLING RATE FOR EACH D¢#) LOOP 
1 - COMMON SAMPLING RATE FOR ALL DCZ) LOOPS 

2 - INDIVIDUAL SAMPLING RATE FOR EACH PID LOOP 
3 - COMMON SAMPLING RATE FOR ALL PID LOOPS 


PROGRAM START ADDRESS: @@090@2 
PROGRAM RESTART ADDRESS: 990063 


TELETYPE OPERATION CODES 

CTRL/T - ENTRY OF SAMPLING RATE 

CTRL/D - ENTRY OF DCZ) PARAMETERS 

CTRL/P - ENTRY OF PID PARAMETERS 

CTRL/C - ENTRY AND EXIT OF CALCULATOR MODE 


REAL TIME CLOCK FREQUENCIES 

Pree POM Haw 

2 - 100 HZ. 

3 - 1000 HZ. 

SAMPLING RATE IS N*(€1/F)5 WHERE N=INTEGER 


A/Ds D/A - ANALOG COMPUTER LINKS 
COMPUTER @ - USES A/D» D/A CHANNEL @ 
COMPUTER N - USES A/Ds D/A CHANNEL N 


PROGRAM SERVICE SEQUENCES 
COMPUTER @, COMPUTER 1; coe »COMPUTER N 


CALCULATOR MODE PERFORMS: 

A+Bs A-Bs A*Bs OR A/B 

ENTER: NUMBER As DESIRED OPERATIONs NUMBER Bs CRe 
APPROXIMATE RANGE OF MAGNITUDES IS 7-2*10**75 TO 8-7*10**-78 


ENTER CONTROL AND CALCULATOR PARAMETERS 
IN ENGINEERING NOTATION. CAUTION: AN INDEFINITE 
STRING OF CHARACTERS WILL BE ACCEPTED. 


THERE ARE FOUR D/A CHANNELSs D/A OVERFLOWS 
WILL BE SET TOige7= 42 VOLES 

MAXIMUM NUMBER OF CONTROL LOOPS: 8 

THIS CAN BE EXPANDED BY INCREASING THE 
PARAMETER STORAGE AREA 


A TIMEOUT WILL OCCUR IF ALL CONTROL LOOPS 
ARE INCOMPLETE IN REQUIRED TIME 


RESPOND TO Ns: WITH A DIGIT REPRESENTING DESIRED OPERATION 


Figure B.1 Start-up Operating Message 
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Figure B.2 Hybrid Computer System 
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